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ABSTRACT 

We use a new method to model fluctuations of the Lyman- Werner (LW) and Lyman- 
a radiation backgrounds at high redshift. At these early epochs the backgrounds are 
symptoms of a universe newly lit with its first stars. LW photons (11.5-13.6 eV) are of 
particular interest because they dissociate molecular hydrogen, the primary coolant in 
the first minihalos. By using a variation of the halo model, we efficiently generate power 
spectra for any choice of radiation background. We find that the LW power spectrum 
typically traces the matter power spectrum at large scales but turns over at the scale 
corresponding to the effective 'horizon' of LW photons (~ 100 comoving Mpc), unless 
the sources are extremely rare. The series of horizons that characterize the Lyman-a 
flux profile shape the fluctuations of that background in a similar fashion, though those 
imprints are washed out once one considers fluctuations in the brightness temperature 
of the 21-cm signal. The Lyman-a background strongly affects the redshifted 21-cm 
signal at just about the time the LW background begins to dissociate H2, so measuring 
that background's properties will reveal important information about the transition 
from early Population III stars to more normal stars. Around this time we flnd that 
fluctuations in the LW background are weak; the fractional standard deviation is less 
than ^ 0.5 on scales > 10 cMpc, only rising to be of order unity on scales < 1 cMpc. 
This should not lead to substantial spatial fluctuations in H2 content, except at the 
earliest times. Even then, most halos form far from other sources, so the transition 
from star formation in low-mass to high-mass halos is rather homogeneous across the 
universe. 

Key words: cosmology: theory - first stars - galaxies: haloes - galaxies: high-redshift 
- stars: Population III 



1 INTRODUCTION 

An important aspect of the cold dark matter (CDM) uni- 
verse is that density fluctuations exist on small scales. These 
small scale perturbations are superimposed on larger scale 
perturbations; the density reaches its highest value over 
the smallest region. Consequently, structure forms via hi- 
erarchical buildup. An initially smooth density distribution 
eventually morphs into a web of sheets and filaments. It is 
the overdense junctions of these filaments that we call dark 
matter halos. Further structure development takes place in- 
side these halos, commencing with the first (Population III) 
stars. 

Population III (Pop III) stars illuminated our dark 
universe in its cold youth and from them developed the 
complex environment we live in today. According to hi- 
erarchical structure formation, these stars formed out of 
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metal-free H/H e gas contained in mini halos at redshifts 
z ~ 20 - 30 (|Couchman fc ReesI 1 19861 ). The minihalos, 
with masses around 10^ M©, have virial temperatures less 
than lO "* K - below the t hreshold for atomic hydrogen 
cooling l|Oh fc Haiman|[2"002l '). Consequ ently, the halos have 



to rely primarily on H2 for cooling |Haiman et al 



iTegmark et~aLlll997l : I Abel et al.ll2002l : iBromm et ~ 



.19961 : 
2OO2I ). 



This cooling takes place via collisional excitation (mainly be- 
tween H2 molecules and energetic H atoms) and subsequent 
radiative decay of the rotational transitions of H2 . 

In the classical view of Pop III star formation, molecular 
cooling produces a single, massive star from cold gas that 
becomes trap ped in the dark matter potential well of one 
minihalo (Bromm et al.l [20091 ). A dense core, or protostar, 
gradually emerges and grows into a massive star by accret- 
ing the surrounding gas. These first stars co uld theoretically 
grow to be several hundred solar masse s (iBromm fc Loebl 
2004); most were probably ~ IOOMq (|Bromm fc LarsonI 
200J). But what would happen if the infalling gas became 
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fragmented? Several studies show that a primordial proto- 
stellar cloud will most likely not violently fragment enough 
for a secondary clum p to compete with the parent clump an d 
form a second star (jAbel et al.ll2002l : lYoshida et "al]|2008l ). 
However very recent simulations suggest that if a gas cloud 
surrounding a protostellar core has an initial degree of an- 
gular momentum, it could collapse into a dense disk, cool, 
and fragmen t, resulting in a binary or even mult iple Pop III 
star system (|Stacv et al.ll2O10l : iTurk et al.ll2009l ') consisting 
of two or more lighter stars as opposed to a single, massive 
star. These new studies could indicate that the formation of 
the first stars could be more complicated and varied than 
previously believed. 

In these primordial star cookers, the ability to form new 
stars is dependent on the abundance of II2. However, as 
the population of these luminous stars grows, the sites of 
star formation are increasingly irradiated by soft UV pho- 
tons (the Lyman- Werner, or LW, bands: 11.2-13.6 eV) from 
existing stars. This LW radiation can photodissociate the 
H2 molecules in the gas through the two-step Solomon pro- 
cess (,Field et al.. ,1966: Stcchcr & Wilham& ,1967, ). 

H2 -f 7 2H, (1) 

in which an H2 molecule hit by a LW photon bumps it 
up to an excited electronic state, H2. A fraction of decays 
from this excited state end up in the vibrational contin- 
uum of the ground state, dissociating the molecule. If the 
LW background becomes strong enough, it can prevent fur- 
ther collapse and consequently stall the further formation 
of primordial ionizing sources by terminating the minihalos' 
primary cooling supply (|Haiman et al.. ,19971 ). As a result, 
the only halos able to cool (via atomic line cooling) and 
form new stars are those with Tvir ^ 10^ K, or masses above 

^ 1O*M0[(1 -f 2:vir)/10]"^/^. Minihalos, with temperatures 
below the threshold for atomic cooling, will not be able to 
collapse past virialization without a sufficient supply of H2 . 

Although these larger halos may still form metal-free 
stars, the thermodynamics of the cooling process is suffi- 
ciently different that we expect the resulting stars to dif- 
fer substantially (especially in their characteristic mass). 
The two populations are sometimes described as Popula- 
tion III.l and Population HI. 2 to emphasize this: both may 
be 'primordial,' but they have very different properties reg- 
ulat ed by the LW (and other) radiation backgrounds (see, 
e.g., iBromm et a"l]|2009l ). 

Most previous calculations of the LW background used 
a homogeneous approximation, in which they as s umed 
a uniform distribution of sources (IHaiman et al ] I2OOOI : 
iRicotti et al.ll2002l : lYoshida et aLlbOOal ). But the highly clus- 
tered, discrete sources responsible for the background radi- 
ation do not generate a uniform background. If these fluc- 
tuations are large enough, the transition from H2 cooling 
to atomic line cooling would be very patchy, potentially al- 
lowing exotic star formation to persist for long periods even 
after the mean background reaches t he threshold value for 
H2 suppression. iDiikstra et al.l (|2008t ) were the first to con- 
sider the inhomogeneous LW background, but only in the 
context of close halo pairs (using a Monte Carlo model) and 
only when the backgr ound was already well above thresh- 
old. lAhn et all (|2009D were the first to consider the inho- 
mogeneous background using a large-scale radiative transfer 
simulation of reionization. 



These photons have other observable efl'ects as well; 
most importantly, as they redshift into the Lyman-a tran- 
sition they couple the excitation temperature of the 21- 
cm transition of hydrogen to the gas kinetic tempera- 
ture via a radiative pu mping mechanism known as th e 
Wouthuysen-Field effect (|Wouthuvsenl 1 19521 : iFieldl Il959l ). 
This rend ers the 21-cm signal visible in emission or absorp- 
tion. Barkana fc Loe bl (|2005a ) showed that the fluctuations 
in this young Lyman-a background produced strong fluctu- 
ations in the 21-cm signal. Conversely, observing these fluc- 
tuations can reveal a wealth of information as to the prop- 
erties of these first luminous sources. We will see that these 
photons begin to affect the 21-cm background at roughly 
the same background intensity at which they suppress H2 
cooling. Thus redshifted 21-cm measurements offer an ex- 
cellent chance to study the transition from Population III.l 
to HI. 2 stars as well as the inhomogeneities in the ultraviolet 
radiation field during the 'cosmic dawn.' 

In this paper, we present a new method with which to 
efficiently calculate the power spectrum of an arbitrary ra- 
diation field for any desired redshift and range in scale (in 
this paper we f ocus on the LW and Lyman-g backgrounds 
specifically; see iMesinger fc Furlanet"tcl |2009| for an earlier 
application specific to the ionizing background). Using the 
halo model to determine the spatial distribution of halos, 
we can build up the radiation background by superimposing 
a flux profile specific to that particular background on each 
halo. This profile effectively replaces the mass density profile 
traditionally used in the halo model to calculate fluctuations 
in the density field. We aim to study the importance of fluc- 
tuations in these backgro unds and complem ent the radia- 
tive transfer simulation of lAhn et all (|2009l ) with our sim- 
ple, analytic model. Our method also takes a very different 
approach to calculating 21-cm fluctuations (due to pertur- 
bations in the Lyman-a radiation field) compared to exist- 
ing w ork (jBarkana fc Loebl l2005al : iPritchard fc Furlanetto! 
I2OO6I ). 

In this first exploration of the radiation background, we 
restrict our attention to the soft-UV background from a rel- 
atively simple model of first galaxy formation. In fact, many 
other physical factors contributed to the transition from 
Population HI to Population II star formation. The most 
obvious is metal enrichment, which also affects the cooling 
and is highly inhomogeneous (see, e.g., iFurlanetto fc Loebl 
l2005h . Also, X-rays emanating from the first sources can 
counteract H2 destruction by increasing the free electron 
fraction and so catalyzin g its formation (jMcDowelll Il96ll : 
iHaiman et al. 1 11996', 2000"). There has been considerable de- 
bate as to which of these backgro unds is more infiu e ntial. For 
our simple model, we wiU follow iMachacek et Zl (l2003ll by 
assuming that the enhancement of the electron density due 
to the X-ray background occurs too slowly to compete with 
photodissociation and so neglect the X-ray background. 

Recently, iTseliakhovich fc Hiratal l|2010l ) pointed out 
that the residual relative velocities of the baryon fiuid and 
underlying dark matter distribution, imprinted during the 
recombination era by the baryons' close coupling to photons 
and now visible as baryon acoustic oscillations, may have im- 
portant implications for star formation in these early, fragile 
halos. These large-scale velocities wi ll suppress the accretion 
of ga s onto small dark mat t er ha los (|Tseliakhovich fc Hiratal 
I2OIOI : iTsehakhovich et al.l I2OI0I ). The actual implications 
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for star formation are as yet unclear; the first sirnula- 
tions show modest effects on the reionization era itself 
(jMaio et al.ll2010l : Istacv et al.ll2010l ). but the effect on the 
earlier ep ochs is important for the LW and Lyman-a back- 
grounds (|Dalal fc PenlbOld ). Because the effects are as yet 
unclear, we will ignore these velocity corrections here, thus 
providing a baseline prediction for comparison with future 
work better incorporating them. 

This paper is organized as follows: in Section [2] we de- 
scribe our method for calculating the power spectrum of the 
LW and Lyman-a radiation background fluctuations using 
the halo model. In Sections [3] and [4] we calculate the flux 
proflles for the LW and Lyman-a backgrounds, respectively, 
and also present our results. We summarize our results and 
conclude in Section [S] We adopt a background cosmology 
(fio, Oa, ^b, h, as, n) = (0.26, 0.74, 0.044, 0.7 4, 0.8, 0.95) con- 
sisten t with the most recent measurements (|Komatsu et al.l 

[ml). 



2 METHOD 

We are interested in modeling the power spectrum of fluctu- 
ations in the LW and Lyman-Q radiation backgrounds using 
the halo model. Unlike earlier treatments of the LW back- 
ground, we are speciflcally interested in its large-scale inho- 
mogeneities, complementing the high res olution, large-scale 

(|2009l ) 



N-body radiative transfer simulation of lAhn et al. 



and the small-scale treatment of iDiikstra etaP ll2008l). In 



stead we will expand the model of Mesinger fc Furlanettd 
(|2009t ). who treated the inhomogeneous hydrogen-ionizing 
ultraviolet background using a halo model- like presc r iption . 

The halo model, as described in ICoorav fc Shethll|2002h . 
uses properties of virialized dark matter halos to calculate 
the effects of non-linear gravitational clustering, assuming 
that all mass in the universe is compartmentalized in such 
halos, whose properties can be parameterized purely by their 
mass m. The three ingredients of the model are the (1) halo 
number density, n(m), (2) spatial distribution of the halos, 
and (3) distribution of mass within each halo, or halo den- 
sity proflle, p{r\m), (where r is the distance away from the 
center of a halo with mass m). Typically, a theoretically- 
mo tivated halo mass f unction for (1) (the classic choice be- 
ing |PressI&|^SchecMe3[l974.) allows one to calculate (2). The 
density proflle for (3) ca n be calibrated by nu m erical simula- 
tions, such as the NFW ([Navarro et al ]|l996l ) or lMoore etal] 
l|l999l ) proflles. This model is very powerful in that it can 
efflciently determine the power and many other useful prop- 
erties for any density fleld at an arbitrary epoch and scale. 

Since we wish to quantify the radiation background 
rather than the mass density fleld, we simply replace the 
'halo density proflle' with the proflle of the radiation fleld 
around each halo. This flux proflle, pia,d(r\m), depends on 
the radiation background under consideration and will be 
discussed later. For a spherically symmetric proflle, the nor- 
malized Fourier transform, u{k\m), can be written as: 



u{k\m) 



dr47rr^ [sin(A:r ) / (fcr )] p^ad {r\m) 



(2) 



where Vc is the cutoff distance at which an observer can no 
longer see the radiation emanating from the source. For the 
case of the LW radiation, for example (described more fully 



in section [3.1^ . this cutoff distance, or horizon, is given by 
''LW 100 comoving Mpc (cMpc). 

As described above, we are interested in modeling the 
fluctuations in a variety of radiation backgrounds (in this 
paper, the LW and Lyman-a backgrounds). We use the 
power spectrum, P{k), or the dimensionless quantity A(fc) = 
P{k)/2Tv^ to quantify these fluctuations. Following the 
halo model, we write the power as a sum of two terms: the 
first term, P^^{k), describes the case for which radiation at 
two points comes from the same sourceQ while the second 
term, P^^{k), describes the case for which the two points 
are illuminated by two sources: 



P(fc) 



P"^(fc) 



P^'^(fc) = 



P^^{k) + P^^{k), where 



dmn{m) 



dmn{r. 



P/coll 



P/c< 



|u(fc|m) 



u{k\m)h{m) 



(3) 

(4) 



Here, p is the matter density, /coii is the collapse fraction 
(fraction of mass in the universe contained in galaxies, or 
collapsed in a halo), b(m) denotes the halo bias (describin g 
how strongly clustered the halos are; iMo fc Whit'3ll996l ). 
and lastly P''°(fc) is the linear power spectrum. Here we have 
approximated the halo-halo power spectrum, for two halos 
with mass m\ and m2, as ti(mi)6(m2)P''"(A;), which requires 
that the halo fluctuations remain linear on the appropriate 
scales (i.e., those on which P^^ dominates). While the den- 
sity fluctuations themselves are very weak at the redshifts of 
interest to us, the halos are also highly biased, so nonlinear 
correcti ons will be important on sufficiently small scales. We 
use the lEisenstein fc Hul (|l999l ) fit to the transfer function 
to calculate P''"(fc) and the Sheth-Tormen mass function 
and collapse fraction jSheth fc Tormeiill 19991 ). 

To model the density power spectrum one must include 
the entire halo population, over all masses. However, we are 
interested in the total radiation field and so should not in- 
clude the low-mass halos unable to host stars. We assume 
that only halos more massive than a cutoff mass, Mmin, host 
stars and so contribute to the radiation background. To mo- 
tivate our choices for Mmin, we first consider the 'filter mass,' 
-A^flitcr, the characteristic scale over which baryonic pertur- 
bations are smoothed in linear perturbation theor y or the 
minimum mass of a halo to a ccrete baryons (jGnedin fc Hull 
1 19981 : iNaoz fc Barkanall2007h . as a lower limit (~ IQ^Mq 
in our redshift regime). In linear theory, the relative force 
balance between gravity and pressure can be characterized 
by the Jeans mass, Mj; the corresponding Jeans scale is 
the minimum scale on which a small perturbation will grow 
due to gravity. M,j depends on the instantaneous value of 
the sound speed of the gas, consequently overestimating 
the c haracteristic m ass scale by up to an order of magni- 
tude l|Gnedinl 12000"). In contrast, Mflner, which takes into 
account the full thermal history of the gas, is a more accu- 
rate mass scale. 

An upper limit would be the threshold for atomic 



^ In the context of the halo model, our radiation background cal- 
culation is analogous to the dark matter density power spectrum 
in the halo model, not to the (discrete) galaxy power spectrum. 
Thus the 'one- halo' term is very important to our results on small 
scales, even if each halo contains only one galaxy. 
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10 100 
R [cMpc] 

Figure 1. The normalized LW flux proflle, PlwC*"!™)! shown for 
2 = 10, 25 (right and left sets of curves, respectively). Each profile 
terminates at the horizon tlwi indicated by the vertical dashed 
lines. The dotted curves represent the corresponding time depen- 
dent versions of the flux profile using the timescale of typical halo 
growth, U, for Mmin = lO^M©. 



cooling, Tvir ~ 10^ K (|Oh fc HaimanI |2002| ). Since the 
H2 f raction, /h, , i ncrea ses with halo mass (/h2 oc 
T^i^ , iTegmark et al.l llQOTf ) the cutoff mass certainly lies 
somewhere between these two limits. The classical criterion 
that the cooling time be smaller than the dynamical time 
will set the redshift- dependent transition: in the absence 
of a LW background, these successful minihalos probably 
have /h, - 10"" and Mhaio - 10*^ Mr (iHaiman et al.l 1 19961 : 
iTeemark et al.l Il997l : lYoshida et"al] I2OO3I ). However, rather 
than try to model this in detail we will employ a variety of 
selections for Afmin in order to remain most general. 



3 THE LW BACKGROUND 

In this section we will apply the above method to fluctua- 
tions in the LW radiation background, which determines if 
the sterilization of minihalos at high redshift (through the 
photodissociation of H2) was a patchy or homogeneous tran- 
sition. 



3.1 The Flux Profile 

Our LW flux profile (shown in Figure [T] for z = IQ and 25) 
for a halo with mass, m, located at an effective luminosity 
distance, r, from the observer is given by: 



Prad(rlm) 



^ /mod(r) 



where we assume for simplicity that the luminosity of each 
halo scales with its mass. Her e we use t h e pick et fence 
modulation factor, /mod(''), from I Ahn et al.l l|2009l ). This is 
the fraction of LW continuum radiation emitted by a source 
that is received by the observer without redshifting into a 
hydrogen Lyman series resonance line, where it will either 
be absorbed or scattered. An absorbed photon will either 
cascade to the 2p level and produce a Lyman-a photon or 
cascade t o the metastable 2s level and decay by two photon 
emission (jPritchard fc Furlanett"oll2006l ')- either way, the re- 
sulting photon will be below the LW range. On the other 
hand, the scattered photon will be reabsorbed until it, too, 
decays into a low-frequency photon (typically after just a few 
scatterings). So, for a given photon at observed frequency, 

correspond- 



Vohs 



we can define a maximum redshift, 



ing to the maximum distance within which photons from a 
source remain in the LW band without redshifting into the 
closest Lyman line from above (located at frequency Vi): 



l + 2o 



(7) 



With each Lyman line associated with its own Zmax.i 
and the spacing between them decreasing with increas- 
ing Vi, we are left with a transmission spectrum resem- 
bling a poorly fas h ioned picket fence, illustrated in Figure 
2 from lAhn et al.l l|2009l ). The modulation factor, /mod, is 
defined as the fraction of the LW frequency interval, 11.5- 
13.6eV, that lies within the pickets, or that is successfully 
transmitted to the observer: 



/mod — 1 — 



V 2.1 eV 



(8) 



where A!/gap,j is the frequency interval between each picket 
in which there is no transmission. The profile terminates 
at the 'LW horizon,' tlw = 97.39q cMpc, the distance at 
which a photon redshifts across the maximum picket spac- 
ing (between the pickets corresponding to the hjS and Ly7 
lines). The scaling factor, a, is defined as: 



h_ 

07 



027 



-1/2 



21 



-1/2 



(9) 



(6) 



W hile /mod can be calculated numerically, lAh n et al.l 
(|2009l) have devised a fitting formula: 

/mod(7-) = 1.7exp [-(rcMpc/116.29a)° ''*] -0.7 

if rcMpc/a < 97.39 and zero otherwise, where rcMpc is the 
distance to the source in cMpc. We have successfully repro - 
duced /mod using the method described bv lAhn et al.l (|2009l ) 
and have confirmed that the fitting formula is accurate to 
within 2 per cent error of the true numerical values. 



3.1.1 The Light Cone 

There is one difficulty with the halo model as usually con- 
structed for our problem: it does not allow the sources to 



^ A more complex relationship is likely, but any such relationship 
can be bracketed by our different choices of A/min- 
^ We implicitly assume a fiat photon spectrum within the LW 
range here, which is a reasonable approximation over this short 
frequency interval. 
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evolve over time. As usually constructed, the halo model 
takes the properties of each halo at a particular instant. 
This is not actually appropriate for our application, where 
the time delay from the finite speed of light implies that 
many sources will only be visible to a given point as they 
were long in the past, when their luminosity may have dif- 
fered from the present value. For example, the light travel 
time across the LW horizon is 38.6 Myr at z = 10 and 11.1 
Myr at z = 25. These values are a full ~ 5% of the Hub- 
ble time at those epochs, so a fraction of the visible sources 
would appear much dimmer than the above model would 
suggest. We next estimate how much we would expect the 
inclusion of such a time dependence to alter our results. 

We can crudely account for halo growth by attaching a 
damping factor to the flux profile: 

pMr\n.) = L{m)^ Lim)h^^g^, (10) 

where — ct*(l + z) in cMpc and corresponds to the 
typical timescale for halo growth, assuming that at these 
high redshifts the halos grow exponentially fast so that the 
luminosity of a halo L{m) oc exp*''** . The growth timescale 
we define as: = a{z)tii (where a is some proportionality 
factor that evolves over time and tn is the Hubble time). 

To estimate a, consider a population of identical halos 
with mass m, at some redshift z, in which the halos are 
conserved; no new halos are created and none are destroyed. 
In this simple case, the collapse fraction is given by: 

. 111^. (11) 
P 

Note how the halo mass is now a function of time, t. Taking 
the time derivative of this expression leads to: 

1 rf/coii _ 2_dm _ J_ 
/coii(t) dt ^ m dt U' ^ ' 

Thus for example, a{z = 10) = 0.283,0.167 for M^in = 
1O'*M0,1O*M0 respectively and a(z = 25) = 0.062,0.034 
for the same choices of M^i^. 

Although this simple model overestimates the growth 
rate of individual sources (because in reality much of the in- 
crease in collapse fraction is driven by new halos passing the 
relevant mass threshold), it provides a simple conservative 
parameterization of the effects of growth. 

Adding this rapid source evolution effectively damps 
the source profile whenever r* "^tlw - these effects are il- 
lustrated in Figure [1] For example, for Mmin ~ 10* Mq, 
r*(2 = 10) = 401.6 cMpc and r*(z = 25) = 53.7 cMpc 
while tlw = 129.7 cMpc and 84.4 cMpc for those respective 
redshifts. We include this crude model for the light cone ef- 
fect below but also point out where it modifies our results 
substantially. 



3.2 The Threshold Intensity 

In order to determine when fluctuations are most important, 
we next compute the evolution of the LW intensity - and 
hence the point at which H2 cooling is suppressed - in some 
simple models of structure formation. This section is not 
meant to provide a detailed model of star formation, but it 
should provide some context for the fluctuations we will later 




10 20 30 40 50 

redshift, z 



Figure 2. The mean intensities of the LW background (solid 
curves), Jlw,21, and the Lyman-a background (dashed curves), 
Ja,2i, for iW^in = Mfiiter, W^Mq, IO^Mq, and 10«Mq (right to 
left). Horizontal dotted lines indicate the upper and lower limits 
for the expected threshold value of Jlw 21 • The dot-dashed line 
indicates the critical intensity for the Lyman-o background, 21 
(see 

examine. We can estimate the mean LW intensity, Jlw (2:), 
with the following: 

Jlw{z) = — J ;^7^^(^')/mod(2' - Z), (13) 

where /mod(^' — z) is part of the LW flux proflle (described 
more fully in >]3.ip and the mean emissivity, e{z'), is given 
by: 

= f.nl^^fcon{z)e,, (14) 

with /* being the star forming efficiency (fraction of baryons 
that actually form stars), which we take to be 10% as a 
fiducial value, and is the mean baryon number density. 
We can approximate the spectral distribution function (de- 
fined as the number of photons per frequency 1/ emitted 
per baryon), etiv), as its mean value et over the LW range 
(11.2-13.6 eV) for simplicity since we are looking at such a 
small range in frequency. We normalize to produce 

4800 photons per baryon between Lyman-a and the Ly- 
man limit for very massive Population III.l (zero-metallicity, 
M > IOOM0) stars (Barkana & Locb 2005a). Note that the 
light cone effect is inherent in this expression due to the 
redshift dependence of the mean emissivity. 

Our results are summarized in Figure (2] The LW inten- 
sity is calculated in units of Jlw, 21 = Jlw/(10~^^ erg s~^ 
cm~^ Hz~^ sr~^). The solid curves represent the very mas- 
sive Population III stars. Replacing the emissivity with that 
of a Population II star (with metallicity equal to 1/20 the 
solar value) p roducing 9690 photons per baryon in our fre- 
quency range (jBarkana fc Loebll2005a^ boosts the intensities 
by a factor of ~ 2; the background reaches threshold earlier. 
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According to iHaiman et al] (|200Cll ). background intensities 
of Jlw,2i ~ 10~^ - 1 are needed to suppress H2 cooling 
in all minihalos over a range of redshift from 2; ~ 10 — 50. 
The lower value describes H2 suppression in halos near the 
Tvir < lO'^ '' K limit; below these low temperatures H2 cool- 
ing is inefficient even in the absence of any photodissociating 
background, so no stars will form. Since /hj increases with 
Tvir, it will be more difficult to terminate a more massive 
halo's larger cooling supply. Thus, as the halo population 
evolves, becoming more numerous and more massive with 
time, the threshold intensity must increase, self-regulating 
star formation by shifting the minimum mass to higher val- 
ues. Once halos reach Tvir > lO'' * K the value of the back- 
ground intensity is once again irrelevant since these large 
halos are able to cool via atomic line cooling and no longer 
rely on their fragile H2 supply. In further calculations, we 
will take Jlw,21 =0.1 as our fiducial threshold intensity. 

It is evident from Figure [5] that our models reach this 
threshold between redshifts 2; ~ 15 — 35. Of course, these 
models are extremely naive and ignore a host of compli- 
cations (such as the evolving star formation efhciency and 
cooling threshold, as well as other feedback mechanisms). 
But they suffice to illustrate approximately when a given 
model reaches the H2 photodissociation threshold, where the 
fluctuations which we will study are particularly interesting. 
Note that, because structure formation itself proceeds expo- 
nentially fast at high redshifts, uncertainties in the star for- 
mation parameters themselves are relatively unimportant. 
In any case, Jlw,2i oc so it is easy to read off the 

appropriate intensity for such a model. 

We emphasize that the fluctuations will be most impor- 
tant near the threshold, because that is when the transition 
in cooling modes actually occurs. If fluctuations are small, 
the transition would occur uniformly over the entire Uni- 
verse. If not, the Universe could contain isolated, sparsely 
populated patches in which H2 cooling remains possible. If 
a minihalo inhabits one of these 'safe' patches, it could con- 
tinue to form very massive stars via H2 cooling even after 
the mean intensity reaches threshold. On the other hand, 
even well before an average IGM point reaches threshold, 
regions near existing sources will be well above it, and this 
could strongly affect the highly-clustered early sources. We 
will examine this phase in i)3.4l 

3.3 Power Spectrum 

Results are depicted in Figures [3] and U In the former, we 
simultaneously vary Mmin and redshift so that Jlw,21 = 0.1 
is fixed, while in the latter we follow a single star formation 
model over redshift (varying Llw,2i). 

The normalized scenarios displayed in Figure [3] include: 
Mmin = 1O^M0, z = 30.15 (bottom panel) and Mmin = 
IO^Mq, z = 20.5 (top panel). The short-dashed curves rep- 
resent the light cone versions for both scenarios, using equa- 
tion (|10p for the flux profile. It is evident that the inclusion 
of the light cone effect preserves the shape of the power but 
modestly boosts the amplitude (by a factor of ~ 2). We 
have also separately displayed the 1-halo and 2-halo terms 
(bottom-most, solid lines) for the Mmin = 10^ M© scenario 
so as to gain a sense of when these terms are dominant and 
to show how they work in tandem to determine the shape 
of A2(fc). 
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Figure 3. Power spectrum of the LW background, A^(fe) for 
scenarios normalized to reach Jlw,21 = 0.1. Scenarios include: 
Mmin = W^Mq, z = 30.15 (bottom panel) and Afmin = 10** Mq, 
z = 20.5 (top panel). The light cone versions for both scenar- 
ios are also shown (short-dashed curves) while the original ver- 
sions are the solid curves. For the Mmin = W^Mq curve we have 
shown not only the total spectrum but also the one-halo (la- 
beled 'Ih') and two-halo (labeled '2h') contributions for visual aid 
(see equation[5ll. The two, topmost, long-dashed curves represent 
the light cone scenarios with the duty cycle taken into account. 
The Afmin = 10^ Mq model includes two different values for this: 
/duty = 0.02 (topmost curve, labeled 0.02) and /duty = 0.06 (sec- 
ond curve, labeled 0.06). The Mmin = 10* Mq case also includes 
two curves for /duty = 0.01 and 0.04. 



It is important to note that the above prescriptions as- 
sume that all halos above the mass threshold, Mmin, form 
stars continuously. Of course, these stars have finite life- 
times, and in the classical Pop III scenario in which each 
halo undergoes only a short burst of star formation, not all 
of these stars are going to be 'turned on' when we take a 
snapshot of the fiuctuations at a particular point in time. We 
can account for this simply by incorporating a duty cycle, 
/duty, into our calculation. This addition exclusively affects 
the one-halo term in equation [S] since both n{m) and the 
effective /coii (which in this model gives the fraction of halos 
hosting active sources) are altered by a factor of /duty, the 
two-halo term remains unchanged. 

The two topmost, long-dashed curves in both panels of 
Figure |3] represent the Mmin ~ 10* Mq and 10^ M0 scenarios 
including the light cone effect using two different values of 
/duty . A reasonable estimate for /duty would be the ratio be- 
tween the average lifetime of a Pop III star, r, and the Hub- 
ble time, tn- The lifetime of these massive stars is believed 
to be a few million ye ars (Myrs) (fBarkana fc Loeb .20oTI : 
iBromm fc Larsonll2004 ). though that remains to be directly 
measured. The topmost curve in both cases assumes an aver- 
age lifetime of r ~ 3 Myrs (/duty = 0.02 for Mmin = 10*^ M© 
and 0.01 for Mmin = 10* Mq) while the second curve assumes 
T ~ 10 Myrs (/duty = 0.06 and 0.04 for Mmin = W'^Mq and 
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Figure 4. Power spectrum of the LW background, A^(fc), for our 
A/min = 10* Mq scenario, which reaches threshold Jlw,21 = 0.1 
at 2 = 20.5. The curves show the power at different redshifts in the 
same model: Jisw,2l ~ 0.001, 0.01, 0.1, 1, and 10 (corresponding 
to z = 25.8, 23.35, 20.5, 17.0, and 11.05 from top to bottom). We 
display the time dependent versions (dashed curves) for the first 
and last of these. 

10* Mq respectively). This greatly increases the importance 
of the one-halo term and so boosts the fluctuations on scales 
below the LW horizon. However, note that the mean back- 
ground intensity also falls by a factor of /duty , so these strong 
fluctuations occur well before threshold is reached. 

The most striking feature common to all power spectra 
is the first turnover located at kuw ~ 0.06 cMpc~^. This 
is a strong signature of the LW flux profile, which termi- 
nates at tlw ~ 100 cMpc (fcLW = 27r/rLw ~ 0.06 cMpc~^). 
Power is smallest for the largest scales and then steadily 
increases until it reaches ^lw- In this regime, regions are 
far outside the LW horizon of each source and so sample 
independent patches in the radiation field. The total power 
is therefore simply proportional to the matter power spec- 
trum multiplied by a mean bias factor (squared). However, 
at /clw, A'^ turns over and begins to fall. This is because such 
scales sample the variations within tlw; if the two points see 
the same halo populations, their radiation amplitudes will 
vary together and the fiuctuations decrease. On the smallest 
scales the power turns up and increases monotonically. This 
indicates where the P^^{k) term becomes dominant, which 
occurs on larger scales (smaller k) for increasing choice of 
Mmin because the sources become more rare. 

Note how the signature shifts to slightly higher k in the 
light cone versions, because the damping scale r* < tlw- 
The turnover is also smoothed out as Mmin increases and 
the one-halo component begins to dominate at larger scales. 
This signature is further smoothed by accounting for the 
duty cycle. The shorter the stellar lifetime (and the smaller 
the duty cycle), the more amplified the one-halo term will 
be relative to the unchanged two-halo contribution. 

Figure [4] shows the Mmin ~ 1O*M0 scenario at several 



different redshifts. Curves are labeled according to their nor- 
malized Jlw,21 values. The central curve (blue in the online 
version) corresponds to the J21 ~ 0.1 normalized version 
(z=20.5). Apparent in Figure [4] is the washing out of the 
turnover at increasing z; at high z halos are more rare and 
the LW background patchier - consequently the 1-halo term 
begins to dominate earlier on scales k < kjjw, thus smooth- 
ing out the key signature. 

One important caveat for our model is the assumption 
of linear bias when computing the 2-halo term in the power 
spectrum. For example, consider the Mmin = 10* Mq model, 
which reaches threshold at z ~ 20. At that time, such a halo 
has b ~ 10. Thus, even though the rms density fluctuation on 
~ 5 Mpc scales is ~ 0.04, the halo fiuctuations are ~ &cr ~ 
0.4, where nonlinear effects are becoming important. The 
steep intensity profiles around these sources make cluster- 
ing somewha t more important for the radiat ion background, 
as found by iMesinger fc Furlanettd (|2009l ). and probably 
enhance the fiuctuations on moderately small scale s by a 
factor of a few. For example, iMesinger fc Furlanettd (|2009l ) 
found from semi-numeric simulations that nonlinear cluster- 
ing tends to smooth out the signature turnover in the power 
spectrum of the ionizing background (where it is due to the 
smaller attenuation length of high- 2: ionizing photons); see 
also the discussion in H3.5l below. 

As discussed in 13.21 we assumed /* = 0.1 here. This 
is likely to be an upper limit, and it could be much smaller 
if, for example, the first star to form in each halo suppresses 
the formation of any others. In this case the radiation field 
would not reach threshold until later, when there are many 
more halos and hence smaller fiuctuations. Our scenarios 
therefore provide upper limits to the fluctuation amplitude 
at threshold. 

Given the unusual shapes of these power spectra, we 
next compute the real-space standard deviation in the in- 
tensity to provide better intuition for the amplitude of these 
fiuctuations. We calculate the fractional standard deviation, 
(t(_R), in the following way: 

a\R) = I dk^P{k)WUk), (15) 

where W^{k) is the 'window function' or smoothing window 
over which we consider varying P{k). We employ a simple 
Gaussian window for computational simplicity: 

WR{k) = e-'='«'/2_ (^g) 

We display (t{R) in Figure [S] for the choice of Afmin = 
10* Mq. From top to bottom, the solid curve is normal- 
ized to 10~^ J21 (z = 30.0), while the others have Jlw,21 = 
10"*, 10-^ 0.01, 0.1, and 1. 

The signature turnover at /clw (see Figures [3] and 2} has 
been lightly imprinted onto the shape of (j{R) in the form of 
a gentle kink at i? ~ 100 cMpc. It is evident that at inten- 
sity levels nearly approaching, at, and beyond the thresh- 
old value, (j{R) is small ( ^ 1) down to very small scales 
(~ 1 cMpc), indicating a fairly uniform background. This 
suggests that it is unlikely for isolated patches still harbor- 
ing H2 to exist and foster star-forming minihalos around the 
threshold. 
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Figure 5. The fractional standard deviation, cr(R), of the LW 
background depicted for Afmin = 10® Mq. From top to bottom: 
the solid curve is normalized to 10~^ J21 (z = 30), the long-dashed 
curve to 10"* J21 (z = 28), the dot-long-dashed curve to 10~^ J21 
(z = 25.8), the short-dashed curve to O.OIJ21 (z = 23.4), the dot- 
short-dashed curve to O.IJ21 {z = 20.5), and the dotted curve to 
I.OJ21 {z = 17.0). We take /* = 0.1 in these scenarios. 



3.4 Fluctuations in the Background at Early 
Phases 

Nevertheless, Figure [5] shows that fluctuations in the back- 
ground are large early on when Jlw is well below threshold; 
a(8 cMpc) ~ 20 for Mmin = W^Mq a.t z = 30. On the 
flip side of asking whether or not scattered H2 driven star 
formation could persist in epochs close to or at threshold, 
these large fluctuations could indicate that even in epochs 
for which the background is substantially below threshold 
there will be patches that are locally at threshold in which 
H2 cooling is suppressed. 

However, in this regime the fluctuations are not gaus- 
sian, so the standard deviation a is not a good representation 
of the importance of the fluctuations. Moreover, because we 
primarily care about the radiation intensity at highly clus- 
tered sites of other halos - where star formation is trying to 
occur - simply taking a pure spat ial average is not nec essar- 
ily the proper approach (see also iDiikstra et al]|2008l ). 

In order to delv e into this new ques t ion, w e follow the 
method presented in iFurlanetto fc Loebl (|2005l ) with which 
they calculated the probability that a collapsing halo forms 
in a region already enriched by galactic winds at high red- 
shift. In contrast, we are interested in calculating the prob- 
ability that a collapsing halo forms in a region with a LW 
background above the dissociation threshold. If the sources 
are very rare, this corresponds to lying within a radius i?thrcs 
of a LW emitting halo. Within this radius, the 'new' halo 
- one that has passed the threshold to form stars - is ir- 
radiated by a local LW intensity above the threshold value 
for suppressing H2 cooling. We can start by calculating the 
fraction of space contained within J?thres of all LW emitting 
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Figure 6. The probability, Qthros (solid curves), that a new halo 
forms within iJthrcs of an existing halo, or within the region irra- 
diated by a local LW intensity above the threshold for H2 suppres- 
sion. Scenarios include Mmin = 10* (rightmost, black curves) 
and Mmiu = IO^Mq (leftmost, red curves). Also displayed is the 
mean LW intensity in units of the threshold level, JLW,2l/'^threai 
using /* = 0.1 and Jthros = O.IJ21. 

halos, (5thrcs('2:), assuming that these regions do not overlap: 

Q^hrcs(2) = dm r]{m)n{'m), (17) 

where 77(m) is the ratio of mass irradiated within i?thros of 
a halo with mass m to that halo's mass: 

m 

For example, for a halo of mass 10* at z = 20 with 
M^iri = 10* M0 and = 0.1, i?thres ~ 3 cMpc and r) ~ 
5.5 X 10*. The corresponding values for 2; = 30 are -Rthrcs ~ 8 
cMpc and ~ 9.7 x 10^. If the flux profile were a pure 1/r^ 
power law, then -Rthrcs oc V'^^ , so rj (x m}^'^ . In reality, the 
modulation factor steepens the fiux profile, so rj is closer to 
flat. 

In the limit in which sources are truly isolated, Quiros 
would be the total filling factor of threshold regions. How- 
ever, as more sources appear, their regions will begin to over- 
lap and - because rj is an increasing function of m - grow 
faster. Because we are only after a crude estimate of this ef- 
fect, we do not worry about overlap here and use Q^hrcs 
as our fiducial estimate. A more sophisticated numerical 
or Monte Carlo model c an easily incorporate this possibil- 
ity l|Diikstra et al.ll2008l ). The model will, of course, break 
down when sources become common enough to sit within 
each other's -Rthrcs (indeed, Qthrcs is not limited to be less 
than unity). Fortunately, in this regime near threshold the 
halo model approach is perfectly adequate. 

Since the newly forming halos are spatially biased and 
preferentially collapse near existing halos, this expression 
is not entirely correct. If we consider two halos, the first a 
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newly formed halo and the second an established LW emit- 
ting halo, then the excess probability that the two halos live 
near each other is quantified by the correlation function, ^gg. 
To linear order, the correlation function can be written as 
^gg = 6newfethrcsC«5, where fcnew = b{Mjnin) aud fcthros are the 
biases of the newly collapsed and LW threshold region re- 
spectively and ^ss is the dark matter correlation function. 
The mean bias of the LW threshold regions surrounding es- 
tablished halos can be written as: 



fethrcs — 



J dmmri{m)b{m)n{m) 
J dmmri{m)n{m) 



(19) 



With this in mind, we can approximate the corrected 
probability that a new halo lives within -Rthrcs of an estab- 
lished halo as: 



;thrcs ^Jthrcs 



1 + 6„ 



s£,ss{Rt 



(20) 



We found that these corrections typically boost Qthrcs by a 
factor of ~ 2 in our scenarios. 

Figure [S] shows Qthrcs for a choice of Afmin = W^Mq 
(rightmost, solid black curves) and Mmin ~ 10* M© (left- 
most, solid red curves) with /* = 0.1. The dashed curves 
in Figure [S] represent the mean intensity relative to thresh- 
old, JLW,2i/Jthres- Qthres iucreases as redshift decreases and 
the LW background builds up and becomes more uniform in 
both scenarios. It is evident from Figure [S] that increasing 
the choice of Mmin delays H2 suppression; increasing Mmin 
from 10* M0 to Mmin ~ 10^ M© decreases Qthres, for ex- 
ample, by a factor of ~ 33 at z ~ 20. By increasing the 
mass threshold one eliminates contributions from a host of 
less massive potential sources, requiring more time for the 
background to strengthen and boost Qthros- Identical calcu- 
lations using /* = 0.01 yield values of Qthrcs that are up to 
factors of ~ 15 smaller. 

Also, note that Qthros increases roughly in proportion 
to Jlw: evidently the mean background provides a good es- 
timate of the volume illuminated by a high intensity of LW 
radiation. However, note that the quantitative similarity of 
this filling factor and Jlw / Jthrcs shown in Fig. [6] is coinci- 
dental and does not occur if we, e.g., change our choice of 
threshold value or /*. 

Once Jlw, 21 reaches threshold, Qthrcs ~ 1.2 (for 
Mmin = 10* M0 and /* = 0.1). At earlier times, the dis- 
crete nature of the sources does substantially increase the 
probability for a new halo to lie within a threshold region 
over and above what one might naively guess from Figure [5] 
for example, when Jlw ~ O.lJthrcs, ~ 15% of the halos lie 
in this regime, even though cr < 1 down to very small scales. 
Nevertheless, we still find that at early times there are rel- 
atively few patches above threshold. 

3.5 Comparison to Other Work 

There has been relatively little work on fluctuations in the 
LW backg round at high red shifts. The most salient compar- 
ison is to lAh n et al.l l |2009l '). who calculated the LW back- 
ground power spectra using a large-scale radiative transfer 
simulation with size Rhox ~ 50 cMpc. They resolved halos 
down to 10* M0 and implemented a simple two-population 
model for galaxies, in which halos with M < 10^ Mq had 
large ionizing efficiencies and larger halos had more modest 
efficiencies. By following the radiative transfer of ionizing 



photons, they also included the suppression of galaxy for- 
mation in halos with M < 10^ Mq following reionization, so 
this higher-efficiency population gradually disappeared. 

Although this simulation box should be sufficiently 
large to include a fair sam ple of the h alo population^ it is 
still much smaller than n.w. lAhn et al.l (2009) therefore used 
a periodic tiling in order to fill out the LW horizon. Unfor- 
tunately, this means that they cannot measure the turnover 
at kjjw, nor the regime in which the two-halo component 
dominates and approaches the straightforward limit b'^P''", 
so interpreting their results is somewhat difficult. 

We are unable to implement the self-regulated reioniza- 
tion model used bv lAhn et al.l (|2009h in our simpler analytic 
model, but we have nevertheless made several test calcula- 
tions to compare our results, choosing a reasonable mini- 
mum mass to match the irs at the different redshifts. Fortu- 
nately, lAhn et all (|2009l ') find that their simulation reaches 
threshold at 2: « 16, when the ionized fraction is less than a 
percent - thus we do not expect the self-regulation to be im- 
portant in the regime of most interest. On the other hand, 
their box has only one radiation source a.t z = 19, so the 
pre-threshold regime probably suffers from their finite box 
size even more than expected. 

Comparing to their Figure 12, we find less power, by an 
order of magnitude or so, in the k ~ 0.3-20 Mpc~^ regime 
probed by their simulations, from z ~ 16-8. (They pro- 
vide only an upper limit at higher redshifts because of the 
many fewer sources in the box.) However, our model does 
help to explain the shape of their power spectrum, which 
(when converted to A^) shows relatively flat power over this 
range, especially at the lower redshifts. This is because it lies 
between the turnover at kuw and the regime in which the 
one-halo term dominates (see, e.g., the lower-redshift curves 
in our Fig. [4]). The flattening becomes more pronounced at 
lower redshifts as the one-halo term decreases in importance, 
thanks to the increased source density. 

Nonlinear clustering is one likely explanation for the 
different amplitudes: the simulations can include the fully 
nonli near clustering of these sources, while our model ignores 
them. iMesinger fc Furlanettol (|2009l ) did indeed find a boost 
of power on comparable scales comparing a halo-model im- 
plementation of radiation fluctuations to semi-numeric sim- 
ulations (in this case, the ionizing background at 2; ~ 6). 
However, they found a much more modest boost (a factor 
~ 2), followed by a steepening toward much smaller scales 
relative to th e halo mod e l pred iction. Reconciling our results 
with those of lAhn et al] l|2009l ) requires a much larger effect. 
One possible explanation is the amplitude of the sour c e fluc - 
tuations, which is ~ 2 times larger in the Ahn et al, (200S) 



model at 2 = 20 than in the IMesinger fc Furlanetto 



20091 ) 



comparison, because the much higher bias of the higher- 
redshift halos compensates for the smaller fluctuations in 
the density field. 

Another possible explanation is the finite source life- 
times imposed in the numerical simulations, which decreases 
the number of sources in the box and increases the im- 
portance of the one-halo term (|Ahn et al.ll2009l : llliev et al] 



* According to the methods of lBarkana fc Loebl |20o3) , the miss- 
ing large-scale modes only suppress the halo population by a few 
percent. 
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albeit in a non-uniform manner within the numerical 
simulation. 

In any case, however, both the simulations and analytic 
models agree that, before the threshold is reached, fluctua- 
tions are relatively unimportant. The large-scale uniformity 
of the background seems robust, in the absence of large- 
scale modulation to the source population itself (as may be 
provide d by relative velocities betwe en baryons and dark 
matter; iTseUakhovich fc Hiratall2010l '). 

In principle, we can also compa re our model with th e 
detailed Monte Carlo simulations of iDiikstra et al] (|2008l ). 
However, they focus exclusively on times far beyond thresh- 
old and very close halo pairs, where our crude approximation 
no longer applies. 



4 THE LYMAN-a BACKGROUND 

The Lyman-a background impr i nts flu ctuations onto the 
21-cm signal (|Barkana fc Loebl 2005al) by way of the 
Wouthuysen-Field effect (|Wouthuvsenlll952l : iFieldlllQSSD in 
which the two hyperfine states of neutral hydrogen are mixed 
via the absorption and reemission of a Lyman-a photon. 
Once the first sources in the universe turn on and amalga- 
mate into a Lyman-a background, this effect drives the spin 
temperature, Ts, to the gas temperature, Tk, resulting in 
a nonzero brightness temperature relative to the CMB, Tj,, 
and al lowing the 21-cm line t o become visible. We can write 
Tb as iFurlanetto et al.ll2006l ): 



Ts 



H{z)/{l + z) 



mK, 



(21) 

where xui is the neutral fraction, (1 + 5) is the fractional 
overdensity of baryons, Tj{z) = 2.73(1 + z) K is the bright- 
ness temperature of the CMB, and dwn/drn is the gradient 
of the proper velocity along the line of sight. We can relate 
T^/Ts to T^/Tk with jFurlanetto et al.ll2006l '): 



1 



Ts 



1 + Xc + Xo 



T-, 



Tk 



(22) 



where Xa and the coUisional and Lyman-a scattering 

coupling coefficients. 

Observing the 21-cm signal could provide us a window 
with which to investigate properties of the exotic sources 
that collectively shaped the nature of this background. This 
means that the Lyman-a background is a directly observable 
effect of (nearly) the same photons that make up the LW 
background. This allows us to measure that feedback process 
directly rather than having to infer it from modeling star 
formation in halos, which is quite difficult. 

4.1 The Flux Profile 

Our construction of the L yman -a flux profile fol- 
lows iPritchard fc Furlanettol l|2006l '). To calculate the 
Lyman-a flux originating from a particular source one must 
consider contributions from all Lyn levels. After absorp- 
tion (and ignoring recombinations directly to the ground 
state, which just regener ate the original photons), a frac- 
tion, /recycle l|Hiratall2006f ) ■ of Lyn photons will be converted 
into Lyman-a photons via cascades through a series of ra- 
diative transitions and will contribute to the total observed 
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Figure 7. The Lyman-a flux profile shown for z = W (black, 
top curve) and 25 (red, bottom curve) with M^\^ = 10* Mq. The 
light cone versions are displayed as dotted curves and the cor- 
responding dashed vertical lines indicate the Lyman-a horizons. 
The z = 10 curve is displaced by a factor of 10 to make it more 
visible. 



Lyman-a flux. The remainder will wind up in the metastable 
2s conflguration and decay via two-photon emission, result- 
ing in no L yman-a photon. We use the v alues of /recycle 
presented in IPritchard fc Furlanettol (|2006[ l. For example, 
/recycle ~ 0.2609, 0.3078 for n = 4, 5 respectively. At large 
n, /recycle asymptotcs to a value of ~ 0.36. However we note 
that quantum selection rules (AL = ±1) forbid a Lyman-/? 
photon from producing a Lyman-a photon. 

Since a source can only be separated from the observer 
by a finite distance before its photons redshift into their 
nearest Lyn transitions, a photon received at redshift z as a 
Lyn photon must have been emitted below redshift Zmax: 



1 + Zmax(n) = (1 -I- 2) 



[l-(n+l)-2] 



(23) 



(l-n-2) ■ 

This imprints a set of horizons on the fiux profile; hori- 
zons become smaller and smaller as you consider higher Lyn 
levels. The series of horizons results in a step-like structure 
of the overall profile. We assume the flux from one halo takes 
the following form: 



Fi 



Lya 



47rr2 



eb,a(j^n)/recycle(n-) 



(24) 



where eb,a(i') is the spectral distribution function (deflned 
as the number of photons per baryon emitted at frequency u 
per unit frequency) described by a power l aw eb{v) oc v"''^ . 
We us e the values for as presented in iBarkana fc Loebl 
(l2005al ) for massive Population HI stars and Population II 
stars. A photon emanating from the source at emission fre- 
quency, v'„, located at redshift z' is absorbed by level n at 
redshift z: 



(1 + ^') 



(25) 
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The sum is ultimately truncated at rimax = 23 to exclude 
levels for which the horiz on lies within the HII r egion of a 
typical (isolated) galaxy (|Barkana fc Loebll2005b^ . 

The Lyman-Q flux profiles for z — 10 (top curve) and 
25 (bottom curve) are displayed in Figure [T] We normalize 
the curves arbitrarily here in order to focus on the shape as 
a function of redshift. The Lyman-a horizon distance, r^ya, 
corresponds to the distance over which a Lyman-/? photon 
would redshift into the Lyman-Q resonance. This is the max- 
imum range a photon can travel and become a Lyman-a 
photon; at z = 10 and 25, TLya ~ 390 and 254 cMpc respec- 
tively. In comparison, the LW horizons for those redshifts 
are ~ 130 and 84 cMpc respectively. We can therefore ex- 
pect smaller fluctuations in the Lyman-a background than 
in the LW background. However, note that the difference is 
not as large as one might otherwise expect because the delay 
from the light travel time already reduces the importance of 
distant sources. 

The dotted curves represent the light cone corrected 
profiles; we treat the light cone effect in the same fashion 
that we amended the LW profile in i]3.1.1l As can be seen 
in Figure the light cone curves begin to diverge from the 
original version on scales 7? > 10 cMpc. The final and largest 
'step' (and more and more of the smaller steps as you look 
at higher redshift) in the profile is effectively beveled as the 
flux begins to prematurely slope downward until it runs into 
the horizon. This difference morphs the shape of the power 
spectrum as discussed further in M.3\ 



4.2 The Mean 21-cm Background 

We next estimate the redshifts for which the Lyman-a back- 
ground fluctuations are important wh e n obs erving the 21- 
cm signal. Following iFurlanetto et al] ()2006l 'l we can write 
the fractional variation of the brightness temperature of the 
21-cm line, ^21, in the following way: 



^21 — PSb + PaSa — 5e„ , 



(26) 



where 5b is the perturbation in the baryonic density, 5a is 
that for the Lyman-a coupling coefficient Xa, and 5g^ is that 
for the line of sight peculiar velocity gradient. The expan- 
sion coefficients, /3i, and their evolution over time determine 
the epochs for which the various perturbations influence the 
fluctuations in Tt- In particular. 



and 



/3 = 1 + 



/3c 



3;tot(l + Xtot) 



a;tot(l -|- 2;tot) 



(27) 



(28) 



Pa is basically the f ractional contribu ti on of the 
Wouthuysen-Field effect IWouthuvsenl Il952l : iFieldl \l95i ) 
to the coupling, where Xtot = Xc + Xa and Xc and Xa 
are the coupling coefficients for collisions and Lyman- 
a scattering. For simplicity in our calculations, and for 
easy comparison to earlier w ork ((Barkana fc Loeb ,200 5a:; 
IPritchard fc Furlanettoll2006l '). we ignore all fluctuations ex- 
cept for those due to density {P5b) and the Lyman-a back- 
ground [PaSa)- We neglect perturbations in the neutral frac- 
tion {Px5x) and the gas kinetic temperature, Tk {Pt5t)- 
The coUisional coupling coefficient was calculated as 



in lFurlanettd l|2006l ). The Lyman-a coupling coefficient can 
be written as 



'J OL 



(29) 



where Sg is a correctio n facto r of order 

unity JChen fc Miral da-Escudel |2004| : iHiratal I2OO6I : 
IPritchard fc Furlanetta I2OO6I ) that we neglect in our 
simple model and Ja is the mean Lyman-a intensity. In a 
similar fashion to equation (|13|) . J a is given by: 



Ja{z) = 



Jrccyclci^^j 



47r H 



(30) 



The critical intensity, j;^_2i = 0.66[(1 -I- z)/2Q] (in 
the units of J21), corresponds to the threshold level of 
J a for which Ts sticks to Tk (|Furlanetto et al.l I2OO6I : 
IChen fc Miralda-Escudj|2004h . How does this threshold in- 
tensity compare to the LW intensity at which H2 cooling 
is suppressed? We have displayed Ja in units of J21 as the 
dot-dashed line in Figure [21 The dashed curves in Figure [2] 
represent the calculated average Lyman-a intensities for sce- 
narios with Mmin = Mflitor(2), IO'^Mq, IO^Mq, and 10* M© 
and = 0.1. Notice how these intensities are larger than 
their LW counterparts (solid curves); the Lyman-a horizon 
distance is a factor of ~ 3 times larger than the LW horizon, 
which not only allows the Lyman-a background to build up 
more quickly but also allows for a more uniform background, 
as discussed in i; 14.11 

Around the time that the LW intensity reaches thresh- 
old for H2 suppression, Ja is also somewhat higher than 
that, and hence very close to J a- As a result, the 21-cm 
background is directly sensitive to the physics of cooling; 
around the time when Ts sticks to Tk numerous minihalos 
are shutting down stellar production as their H2 supplies are 
wiped out. Conveniently, this makes the 21-cm background 
a nearly direct probe of this very interesting epoch in the 
history of galaxy formation. 

We present Pa{z) in Figure [8] for Mmin = 10* and 
/* = 0.1. We find that, for Mmin = 10* Mq, Pa peaks at 
z ~ 22 and is significantly nonzero from 2 ~ 15 — 30; fluctu- 
ations in the Lyman-a background are important over this 
range. For this scenario, the mean LW background reaches 
threshold (Jlw,21 ^ 0.1) by 2 21. 



4.3 Results 

Results for the Lyman-a radiation background power spec- 
tra are displayed in Figure]^ For ease in comparison to our 
LW results, we use the corresponding scenarios from FigO 
normalized to reach Jlw,2i ~ 0.1. Values for Pa in these 
scenarios are 0.89 for Mmin = 10^ M0 at z = 30.15 (bottom 
panel) and 0.83 for Mmin = 10* M© at 2 = 20.5 (top panel). 
The light cone versions are also shown (short-dashed curves; 
original versions are the solid curves). The two, long-dashed 
curves at the top of both panels represent fluctuations for 
these scenarios including the light cone effect and duty cy- 
cle. The topmost curve in both cases assumes an average 
lifetime of r ~ 3 Myrs (/duty = 0.02 for Mmin = IQ^Mq and 
0.01 for Mmin = 10* Mq) while the second curve assumes 
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Figure 8. The Lyman-a coupling coefficient, I3a{z), for M^i^ = 
-'^^filtcr (dasfied line) and 10* (solid line), assuming /* =0.1. 
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Figure 9. Power spectrum of the Lyman-o background, A^(A;), 
for scenarios normalized to reach Jlw 21 =0.1 for ease of compar- 
ison with Fig. [3] Scenarios include: M^nin = W^Mq a,t z = 30.15 
(bottom panel) and Mmin = 10* Mq a.t z = 20.5 (top panel). 
The light cone versions for both scenarios are also shown (short- 
dashed curves) while the original versions are the solid curves. 
The two, topmost, long-dashed curves include the light cone ef- 
fect and duty cycle for the Mj^in = 10* Mq case (/duty = 0-01 and 
0.04; labeled 0.01 and 0.04 respectively) and the M^in = 10^ M© 
case (/duty = 0.02 and 0.06). 



r ~ 10 Myrs (/duty = 0.06 and 0.04 for Mmin = l(f Mq and 
10* M© respectively). 

As discussed in § 14.21 the mean Lyman-a? intensity, Ja , 
is very nearly the critical intensity, , around the time that 
the LW intensity reaches the threshold level for H2 suppres- 
sion (this can be seen in Figure [2)). However, Ja oc /duty: 
stars are 'turned on' for a smaller fraction of the time and 
thus build the radiation background more slowly. Thus the 
duty cycle curves in this Figure are not at the coupling 
threshold. Obviously the fluctuations are boosted on small 
and mid- range scales (see Figure|9]), but this is not surprising 
given that we are no longer probing the threshold epochs. 

Present in the original models are a series of sequen- 
tially damped wiggles, in contrast to the smooth transition 
of the LW power. These result from the (Fourier transform of 
the) discontinuous horizon steps present in the Lyman-a flux 
proflle. Unfortunately, these signature wiggles are smoothed 
out once the hght cone effect is applied to the models. As 
discussed earlier in i]4.1l the light cone effect bevels out the 
horizon steps that give the flux profile its distinctive shape, 
resulting in a more featureless profile and producing a nearly 
featureless power spectrum. 

The power turns over at roughly k^ya ~ 27r/rLya, where 
TLya is the Lyman-a horizon distance (discussed above in 
i]4.ip , except for the light cone versions whose turnovers shift 
to somewhat higher k. In addition, the amplitude of the 
light cone Lyman-a power is roughly a factor of 2 smaller 
than the corresponding amplitudes for the LW background 
for scales of fc ~ 0.1 cMpc~^. This is likely a symptom of 
the larger Lyman-a horizon distances, which are 3rLw 
for these redshifts. Furthermore, the hght cone effect on the 
Lyman-a power is stronger than the corresponding effects 
on the LW fiuctuations. The Lyman-a models corrected for 
halo growth over time are boosted in amplitude by a factor 
of ~ 7, while the LW models receive a boost by a factor of 
~ 1.5. The large Lyman-a horizon allows points to 'see' more 
halos, bolstering the light cone effect on the most distant 
sources. 

4.4 The 21-cm Signal 

Finally, armed with the fluctuations in the radiation back- 
ground and the fluctuations in the baryon density (com- 
puted with the linear power spectrum on these scales), we 
can estimate the 21-cm signal itself. Referring back to equa- 
tion (|26|l . we consider fluctuations in Ti, sourced by per- 
turbations in the matter density, Wouthuysen-Field cou- 
pling (or the Lyman-a flux), and radial velocity gradient 
of the gas. All of these fluctuations are isotropic except 
for the velocity fluctuation, which introduces an anisotropy 
to the power spectrum and can be written as Sg^ (k) = 
—fj-^Si, ([Bharadwai & Ali 2004), where /i is the cosine of the 
angle between the wavenumber k of the Fourier mode and 
the line of sight. This enables the t otal power spectrum for T i, 
to be separated into powers of fi^ (|Barkana fc Loebll2005bl ): 

Pt, (k) = M"Pp4 (k) + (k) + P^O (k). (31) 

The fi^ term, which can be written as l|Barkana fc Loebl 

l2005bl ) 

P^2{k) = 2fi^[pPs{k)+PaPs-c{k)], (32) 
contains contributions of density-induced fluctuations in 
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the Lyman-a flux, where Pg^aik) is the cross-power spec- 
trum for the matter density and Lyman-a radiation back- 
Rrou nd. In our halo mode l, this is very easy to calculate 
fc.f.. lCoorav fc Sheth|[20o3 ): 

Ps~a{k) = Ps-a(k) + Pi'laik), where 



Pslc^ik) 



PtUk) 



dmn{m) 1 — — 

M„,i„ \P./coll 



(33) ^ 



— \us(k\m)\\ua(k\iiM) 
P 



dmn(m) ( — ) us(k\m)b(m) 
VP 



(35) 



dmn{m, . 

M„i„ \P/coll 



Ua(k\m)b{m) 



pi.. 



where it^ is the halo density proflle that describes the dis- 
tribution of mass within each halo. 

Therefore, P^j2 (fc) can easily be used to investigate fluc- 
tuations in the Lyman-Q background at high redshift. We 
display P^2(fc) in Figure [TOl for Mmin = lO^M© at redshifts 
z = 11.05, 17.0, 20.5, 23.35, and 25.8 (from top to bottom) so 
as to correspond to the scenarios presented in Figure U We 
find that fluctuations in Tb increase with decreasing redshift 
and decrease with scale. The increase in amplitude levels 
off once Ja reaches the critical intensity, (which, for the 
scenario in Figure [TOl occurs around z ~ 18). At this point 
the Lyman-a coupling is saturated (a^tot S> 1) and the 21 
cm fluctuations become insensitive to the fluctuations in the 
Lyman-Q background. One can also pick out the 'one-halo' 
term kicking in on small scales for these later epochs (the 
z — 11.05 curve in Figure [TO)) . 

The amplitude and overall shape of our 
po wer spectra are com p arable to those presented 
in iPritchard fc Furlanettd (|2006l ). who used the same 
model to describe the Lyman-a flux but calculated Pjj2 
using a linear transfer function. The amplitudes from 
our m odel also agree with those from iBarkana fc Loebl 
l|2005al ). although the shape of the power differs because 
they neglected the effects of atomic cascading (by assuming 
/recycle = 1) in their calculations. 

Unfortunately, these power spectra do not display any 
sort of tell-tale signature feature such as the distinct LW 
turnover as can be seen in Figures [3] and U or the wiggles 
from the horizon steps, because the density-induced fluc- 
tuations wash them out. This ag rees with previous work; 
recently, IVonlanthen et al] ([20111 ) showed with a radiative 
transfer numerical simulation that the horizon steps present 
in the Lyman-a flux profile left imprints in the differen- 
tial brightness temperature profile just after the first lumi- 
nous sources turned on. However, as time progressed and 
more new sources began contributing to the Lyman-a back- 
ground, the steps were effectively wiped out. 



5 SUMMARY & CONCLUSIONS 

In this paper, we have modeled ffuctuations in the LW and 
Lyman-a radiation backgrounds using a variation of the halo 
model. First, we calculated the LW power spectrum and 
found that the power is characterized by an abrupt cut- 
off at the LW horizon distance, rLw; the power turns over 
at the horizon wavenumber, fcLw, unless the sources are so 
rare that the one-halo term dominates (i.e., correlations are 
determined by the flux profiles of individual sources). 




0.01 - 



Figure 10. Power spectrum of Ti, for the 21-cm transition, 
|T6|[fc3p^2(A:)/27r2]V2 (in mK), for scenarios with Mmin = 
10** Mq for ease of comparison with Fig. |4] Scenarios include: 
z = 11.05, 17.0, 20.5, 23.35, and 25.8 from top to bottom. 



We found that the fiuctuations in the background are 
weak and should not lead to substantial spatial fiuctuations 
in H2 content. Once a population of low-mass halos produces 
enough stars to generate a threshold LW background large 
enough to destroy their own II2 reservoirs used for cooling, 
star formation can only proceed in larger halos with more 
substantial reservoirs. Our model predicts that, by the time 
this threshold is reached, fiuctuations in the intensity field 
will be quite small, so this transition will be rather homoge- 
neous across the entire Universe. 

Though we found fluctuations in the background to be 
small around threshold, on the flip side we also found them 
to be large in those early epochs during which the back- 
ground was approaching threshold. This could indicate the 
presence of patches of collapsing halos that are locally above 
threshold and H2 suppressed. Taking into account the bias 
of sources, we crudely approximated the probability that a 
new halo lives in such a patch and found that even though 
source clustering substantially increases this probability, it 
is still relatively small. Thus, at this time, the majority of 
newly-forming halos can continu e to cool via H2, even in the 
presence of established galaxies (|Diikstra et al.|[2008l ). 

Eventually the first stars will build up the LW back- 
ground by enough to suppress their own star formation; 
the process will terminate when only those halos above the 
atomic cooling threshold can form stars. However, the mode 
of star formation in these halos is very different from the 
minihalos that produce the first stars, so the transition to 
higher-mass halos has i mportant consequence s for the global 
star formation history (|Oh fc HaimanI [20021 ) . This transi- 
tion from Population HI to Population II star formation is 
of course extremely complex, and we have examined only 
a small part of it. The formation of very massive Popula- 
tion HI stars requires two physical conditions: (1) metal- 
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free eas and (2) a reservoir of H2 that allows the gas to 
cool ( Bromm fc LarsonI |2004| ). We have not examined the 
first condition, but the slow speed of galactic winds (com- 
pared to the speed of light) guarantees that metal enrich- 
ment will be very inhomogeneous, and pristine pockets of 
gas could persist until very low redshifts (jScannapieco et eiD 
l2002l : iFurlanetto fc LoeblboOSi ). In contrast, the LW back- 
ground is spatially uniform and so will induce a rapid, homo- 
geneous transformation in the fundamental processes of star 
formation, even when metal enrichment remains inhomoge- 
neous. This will induce a shift from very massive Population 

111.1 stars to less massive - but still primordial - Population 

111. 2 stars that require atomic cooling. 

Next we considered the fluctuations in the Lyman-a 
background in a similar fashion. The Lyman-a flux profile 
imprints a series of wiggles on the shape of the power, corre- 
sponding to the series of horizon steps that characterize the 
profile. Unfortunately, unlike the LW case, these signature 
wiggles are washed out once we account for halo growth over 
time. We found that the amplitude of the Lyman-a power is 
smaller than that for the LW background by a factor of ~ 2 
for scales larger than ~ 15 cMpc. The smaller fluctuations 
are due to the large Lyman-a horizon distance, r^ya ~ 3rLw, 
allowing the halos to 'see' further. 

We used our model for the fluctuations in the 
Lyman-a background to generate power spectra for the 
brightness temperature, Tt, of the 21-cm signal. We 
find our v alues to be in good agreement with previous 
estima tes jPritchard fc Furlanetto! I2OO6I : iBarkana fc Loebl 
l2005al ) that used a very different approach to estimate the 
radiation field fluctuations. We do not see a signature fea- 
ture present in the power in contrast to the distinct LW 
turnover. 

Our relatively simple model, though convenient and 
efficient, has a number of caveats that compromise accu- 
racy. Most notably, we neglect nonlinear effects on the back- 
grounds by relying on a linear approximation for the halo- 
halo correlation function. Thus, our models are only good 
down to the scales for which this linear approximation still 
holds. In addition, we assumed a uniform star forming effi- 
ciency, for all minihalos in our calculations. As the halo 
population grows in size and complexity in the later epochs, 
variations in galactic properties - sourced partly by the sup- 
pression of II2 cooling, but also due to a myriad of other 
factors - will complicate our simple treatment. After com- 
paring our results to the previous simulation from lAhn et al.l 
12009) we find that the shapes can be well matched even with 
the limited dynamic range of the Ahn values (the signature 
turnover is not covered in their range). Although our am- 
plitudes disagree by a factor of ~ 10 at the lower redshifts, 
we both find gentle fluctuations in the background around 
the time it reaches threshold, implying that 'safe' patches 
sheltering isolated sources of H2 are rare by this epoch. 

Furthermore, we have neglected the effects of X-rays in 
our models. If the first luminous sources had hard spectra 
that extended out to X-rays, these X-rays would have far 
reaching effects due to their large mean free paths - they 
could, by catalyzing the formation of new H2, potentially 
countera ct H2 photodissociation by the growing UV back- 
ground ijMcDowelll [l96ll : iHaiman et all Il996l . l2000l ). This 
could stall the transition of Population III to Population II 
stars by allowing minihalos to continue forming new stars. 



altering the make-up of the sources responsible for reioniza- 

tion. 

Finally, iTsehakhovich fc Hiratal (I2OI0I ') argue that 
a new nonlinear effect must be considered in structure 
formation: the supersonic relative velocity between dark 
matter and baryons can suppress the matter power spec- 
trum near the baryonic Jeans scale, altering the abundance 
and clustering properties of the first dark matter halos. 
This effect could be accounted for in a future version of 
our simple model by introduci ng a modulat i on fa ctor to 
the halo mass function. In fact, iDalal fc PenI (I2OI0I ) argue 
that the resulting fiuctuations in the radiation background 
could strongly affect the 21-cm absorption power spectrum. 
However, they assume that the velocity exerts a very 
strong effect on galaxy formation, which may be at odds 
with more detai l ed n umerical simulations (Ma io et al.l 
I2OIOI : IStacv et all I2OI0I ). In any case, such a large scale 
modulation will add more power to the radiation field and 
may have important implications for the homogeneity of 
the LW and Lyman-a backgrounds, since these acoustic 
features appear on comparable scales to the LW horizon 
(~ 100 Mpc). 
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